Soliton excitations as emitted clusters on nuclear surfaces 
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Abstract 



This work reports on calculations of the deformation energy of a nucleus 



> 

in 

. plus an emitted cluster as a soliton moving on its surface. The potential 

barrier against the emission of a soliton is calculated within the macroscopic- 
microscopic method. The outer turning point of the barrier determines lim- 
itations on the geometrical and kinematical parameters for the formation of 
a surface soliton. For large asymmetry, the two-center shell model is used to 
assign a structure to the soliton. Calculations for 248 No with the emission of a 
^ 40 Ca soliton are reported; likewise for 224 Th with the emission of 16 0. Except 

for necked shapes at the very first stages of soliton formation, the greatest 
portion of the deformation path displays rather compact configurations with 
large neck radii. These shape sequences correspond to allowable soliton ve- 
locities. Close to and just beyond the touching point configuration, where the 
shape becomes concave, the width and the velocity of the soliton approach 
zero. The calculations suggest that the emission of a 40 Ca structure is quite 
probable due to a low potential barrier, whereas the emission of an 16 type 
soliton is rather unlikely due to the higher penetration barrier. 
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I. INTRODUCTION 

There are different theoretical models which describe cluster formation and emission 
from nuclei, and most of them use nonlinear partial differential equations. A fundamental 
understanding of nonlinear terms in a nuclear model reveals new phenomena and shapes 
more complicated than linear theory suggests. In this paper we introduce a model for the 
cluster emission, based on soliton solutions of nonlinear equation. 

Soliton structures have been analyzed within the frameworks of hydrodynamics, non- 
linear optics, solid state and plasma physics. Experimental and theoretical results |J sug- 
gest that solitons are non-dispersive, localized waves executing uniform motion that can 
be described by three interrelated parameters: amplitude A, half-width L, and velocity V. 
Furthermore, these structures arise as analytical solutions of non-linear dynamical systems, 
like the Korteweg de Vries (KdV) or Nonlinear Schrodinger (NLS) equations. Soliton the- 
ory has been applied many times in nuclear physics, so far. For exemple, it provided very 
good descriptions of some localized stable surface excitations 0, of spectroscopic factors in 
cluster and a emission [|], and of quasimolecular spectra for a plus heavy nuclei collisions, 
Q. Also, the cubic and quintic NLS equations were used in three-dimensional models for 
cluster emission ||, providing results in very good agreement with TDHF simulations. 

In order to look for the possibility of describing nuclear phenomena such as cluster 
emission || by soliton formation on the nuclear surface, it is necessary to assign a microscopic 
structure to the parent heavy nucleus and the emitted soliton cluster. The microscopic 
substructure further allows one to add shell corrections to the usual macroscopic liquid drop 
energy and thus give a complete descripton of the system, from the initial nucleus with no 
soliton substructure to one with a soliton-like structure on its surface and on out to possible 
cluster emission. 

A straightforward way to accomplish this is to calculate shell effects obtained from the 
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single-particle levels of an asymmetric two-center shell model. One of the centers is placed 
in the middle of a small emitted sphere and the other is the center of the heavy fragment. 
This approach allows a microscopic description of the nuclear evolution from one to two 
independent quantum systems. 

The procedure involves calculating the total potential energy as the sum of the macro- 
scopic energy and shell corrections which is then minimized, which yields a barrier that 
increases as a function of the amplitude of the soliton. Calculations have been performed 
for two possible reactions: 

248 Nq ^ 208p b + 40 Ca 
224 Th ^ 208p b + 16Q 

We select spherical daughters and emitted clusters, in order to fit the asymmetric two-center 
shell model we have constructed, which goes from one sphere (parent) to two necked spheres 
and then to separation. 

In the present model, we describe cluster emission processes by using such soliton-like 
shapes on the nuclear surface of the heavy fragment. For a given cluster geometry, we 
calculate the corresponding soliton parameters (A, L, V) as functions of the separation 
parameter, that is along the static path of the cluster emission process. 

Solitary waves have been shown recently to exist on liquid drops, bubbles, and shells 
@. The non-linear hydrodynamic equations are related with the KdV and mKdV equation 
generating localized patterns ranging from small oscillations to nonlinear ones, including 
solitons. This model has a Hamiltonian structure, and such soliton-like excitations were 
observed experimentally when the shape oscillations of a droplet became nonlinear. It is 
therefore natural to extend that this model to other drop-like systems, from neutron stars 
to hyperdeformed nuclei and fission. 

In the first section we define the deformation space we that we work with and the shapes 
that can be obtained therein. In the second section we give a short description of the 
macroscopic-microscopic model we used in the calculations. Emphazis has been given on 
the asymmetric two-center shell model we constructed in order to approach solitonic shapes. 
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Results for the two reactions given above are presented in the third section. 



II. SPACE OF DEFORMATION 



The problem of describing cluster emission (or nuclear fission) by a convenient 
parametrization of the shape is not new, and it is a decisive factor determining the amount 
of calculations. In the present model we use a new type of parametrization, described by 
soliton shapes. As well as the other parametrizations are in relation with the nuclear system 
and its quantization, the soliton parametrization is related to the KdV and mKdV dynamics 
of the soliton. Moreover, the soliton model takes profite from the high stability in time of 
such shapes. While in general one associates solitons with macroscopic pictures, it is never- 
theless true that the soliton dynamics can be quantized and the KdV (or mKdV) 
equations can be related to quantum systems like those described by Schrodinger equations, 
J5|||], as we noted above. 

Solitons on the surface are described by two asymmetrical spheres smoothly joined one 
to another through a neck region (Fig. 1). There are three independent geometrical param- 
eters which form the space of deformation: the distance R between the centers of the two 
fragments, the emitted small sphere radius R2 and the neck sphere radius R3. The neck 
region is obtained by rolling a sphere of radius R3 around the symmetry axis. Such shapes 
are generated by the following equation written in cylindrical coordinates: 



where the quantities not shown in Fig. 1 are: 23-the position of the center of the neck 
sphere on the symmetry axis, z c \ and z C 2-the positions of the intersection planes of the two 
fragements spheres with the neck sphere, and Z\ and Z2~the positions of the two spherical 
fragment centers. The heavy fragment radius -Ri is obtained from total nuclear volume 
conservation. The soliton solution along the 9 direction 




(1) 
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is characterized by the amplitude A, or the relative amplitude a = A/R\, the half-width L, 
and the angular velocity V. The soliton solutions have a special shape-kinematic depen- 
dence, V ~ A and L ~ 1/y/A, that is, a higher soliton is narrower and travels faster [|TJ. 
This relation can be used to experimentally distinguish solitons from other normal modes of 
excitations (for example by calculating the reciprocal moment of inertia) 0. The amplitude 
of the soliton is related to the two-center shell model by the relation with A = R — R\ + R 2 . 
The halfwidth of the soliton is approximated with 2p(z c2 ), or with the diameter of the cir- 
cular surface within the separation plane between the emitted sphere and the rest of the 
shape 

L 2R > A (3) 

(l + a)(R 1 -R 2 y { ] 

We mention that the cubic and quartic NLS equations are related to the mKdV equation 
by a very simple exponential transform, and actually there is no essential difference between 
the NLS and the mKdV solitons jl]|§. As an example, we notice the connection between 
the KdV equation and the nuclear potential in the Schrodinger equation || , or the relation 
between NLS solitons and coherent states 0. Some applications of the KdV or mKdV- 
solitons are macroscopic, but the soliton solutions can be quantified by standard procedures 

mi- 
ni. MACROSCOPIC ENERGY 

The deformation energy E^f is calculated in a macroscopic-microscopic approach: 

Edef = E macro + 5E s hell + SP. (4) 

The macroscopic part E macro includes the shape-dependent components of the charged liquid 
drop: 
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E m acro — Ec + Ey+E, (5) 

where Eq is the Coulomb energy and Ey+E is the surface or nuclear energy calculated within 



the Yukawa-plus-exponential model [10,11 



The Coulomb energy general form is the double-volume integral: 

E = 1 f f Pe{r 1 )p e {r 2 )d?r 1 d?r 2 . . 

° 2 Jv Jv ri2 

which is split into four parts, two of them being equal to one another |J: 

2 Jvi JVi r 12 Jvi Jv 2 ri2 2 Jv 2 Jv 2 r 12 

where r% 2 = — r 2 \- The first and the last terms represent the self-energies of the two 
fragments, and the middle term is the Coulomb interaction between the two fragments. 
In cylindrical coordinates the three terms are given by: 



Pe r''„ r''„ / 2 zd P 2 



E c = -j- dz dz\ I dip i dip i p — 



Pl- 



[P 2 + 



10 Jz> Jz> Jo r Jo T V 2 9z , 

pp x cos(y? - ipx) + p— sm(<p - ipi) + — 

pi - 2 PPl cos(p - Vl ) + (z- ^a) 2 ]^ 2 , (8) 

where p = p(z, p>) is the nuclear surface equation, and z' and z" are the intersections of the 
surface with Oz axis. 

The general form of the Yukawa-plus-exponential energy is: 

a 2 f f (Tvi \ exp(-ri 2 /a) 3 3 



87r 2 roa 4 Jv Jv V a J r 12 /a 

where r\ 2 = \*i — ra|, a=0.68 fm accounts for the finite range of nuclear forces, and a 2 = 
a s {l — nl 2 ). k is the asymmetry energy constant, and the surface energy constant is a s =21.13 
MeV. In a similar way to the Coulomb part, one obtains three terms: 

87r^roa 4 Jvi M V a / ri 2 /a 



2v/a 2 ia22 r 3 r fr 12 Xexp-r^/a 3 



7rroa 4 Vvi Jv 2 \ a ) r 12 /a 



rf 3 ri / pi _ 2 '/^ rf 3 r 2 (10) 



where a 2 i = a s (l — Klf), I{ = (Ni — Zi)/Ai. For shapes with axial symmetry, each of these 
terms involving a double-folded integration over the nuclear volume can be reduced to a 
three-dimensional integral. 



IV. THE ASYMMETRIC TWO-CENTER SHELL MODEL 



The two-center shell model was developed by the Frankfurt school for symmetric splitting 



I3| and for low asymmetry [ Hfl . Here we present the main steps of the two-center shell model 



we developed for large asymmetry starting from another symmetrical two-center model [TO) 



The general Hamiltonian describing the evolution of the level scheme of the two-center 
shell model is based on two oscillators which split from an initial common oscillator. The 
usual spin-orbit interaction and the l 2 term are constructed as depending on the mass 
asymmetry. Thus, the total Hamiltonian reads: 

H = H osc + V(\ ■ s) + V(l 2 ) (11) 

where H osc is the two-oscillator Hamiltonian, and V(l ■ s) and V(l 2 ) are the spin-orbit and 
the l 2 potentials. 

A. The diagonalization basis 

The oscillator part of the Hamiltonian, H osc , is given in cylindrical coordinates by: 



H„ 



h 



2 



d 2 Id Id 2 d 2 
dp 2 p dp p 2 d(p 2 dz 2 



+ V(p,z), (12) 



2m 

where the asymmetric two-center oscillator potential has the form: 

V{p,z) = -m < (13) 
[ u 2 p2 p 2 + uj 2 2 (z - z 2 f , z>z . 

Here zq is the separation plane coordinate between the two asymmetric systems on the 
symmetry axis Oz. Since we consider the case of two asymmetric spheres, u Pl = u zi = u)\ 
and u P2 = uj Z2 = lo^- 



To obtain an appropiate basis for this system, we consider the intermediate case when 



v P1 — u P2 



oj\. Then the potential is: 



1 | ufp 2 + u{(z + z x f , z < 
uj\p 2 + u 2 (z — z 2 ) 2 , z > 0. 



(14) 



At this point the intermediate Hamiltonian is separable, and one gets the eigenfunctions 
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for the axial degree of freedom and 



exp (imc, 



(15) 



Rn p \m\{p) 



2T(n p + l)aj 



„2 „2 ' 



exp 



(16) 



y\ T(n p + \m\ + 1) 

for describing radial motion where OLi = (mcoi/h), T(x) is the gamma function, and L^(x) 
is the Laguerre polynomial. As one knows, the oscillator energy for oscillations in the plane 
perpendicular on the symmetry axis is: 



Ep^ = huj p (2n p + \m\ + 1) 



(17) 



where u p = lo\. 

Solving the third equation which accounts for oscillations along the symmetry axis, 
we have different solutions for the two regions of the nuclear shape. According to the 
^-dependent potential, 



I | uj\{z + z x f , z < 
V(z) = -mo 

1 u*(z - z 2 ) 2 , z > 



(18) 



where Z\ and z 2 are the centers of the heavy and light spherical fragments, respectively, one 
obtains the differential equations 



d 2 2m E z mlujf(z + zi) 2 
d? + ~h 2 



Tr 



d 2 2m E z m\\o\{z — z 2 ) 2 

d? + ~h 2 r? 



Z Ul {z) = 0,z<0 
Z U2 (z) = 0, z> . 



(19) 
(20) 



s 



At this point it is important to mention that the z = plane is the intersection plane between 
the two systems with u pi = u p2 , whereas the "real" intersection between the asymmetric 
spherical systems is at z = z . The solution for the z-dependent Hamilton equation will be: 

gf(z+zi) 2 ' 



Z v (z) = 



C U1 exp 
C U2 exp 



g|(z-z 2 ) 2 



H Vl [-ai(z + zi)] , z < 
H U2 [a 2 (z - z 2 )\ , z > 



(21) 



where C V1 and C U2 are normalization constants, and H v (z) are the Hermite functions. 

As can be seen from these results, four quantities need to be determined: the two quan- 
tum numbers v 1 and v 2 and the two normalization constants C Ul and C U2 . These quantities 
can be calculated from a system of four equations. From the normalization condition 



/oo 
1^)1^ = 1, 
-oo 



from the continuity of the z-wave function and its derivative at z=0 



Z Ul (z = 0) = Z U2 (z = 0), 



(22) 



(23) 



Z' l/i (z = 0) = Z' U2 (z = 0), 



and from the energy matching condidtion along the 2 axis 



hui(vi + 0.5) = hoo 2 {v2 + 0.5). 



(24) 



(25) 



From these, a basis for diagonalization of the potential differences to obtain the real energy 
values can be calculated. 



B. The asymmetric oscillator system 

Once we have total wave functions, we have to determine differences between the diagonal 
Hamiltonian and the real one. First, the oscillator Hamiltonian has to provide the initial 
oscillator potential when there is only one heavy sphere (starting point). For this initial 
configuration the difference that needs to be diagonalized is 
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AV sphere (z) = < 



±m [ujl(z + zx) 2 - ujI(z - z 2 ) 2 } , z>0 
, z < 0. 



(26) 



AV(z) 



For the next stages of deformation, the difference between the z-dependent oscillator poten- 
tials that needs to be diagonalized is: 

z < 

lm [oul(z + zxf - w\(z - z 2 ) 2 , < z < z ( 27 ) 

, Z > Zq. 

As for the difference in the p-dependent oscillator potential, this only exists for intersecting 
spheres and is given by: 



AV(p) 







, Z < Zq 



\m (ujl - ujl)p 2 , z > z 



(28) 



or, if written as an operator, the quantity to be diagonalized is given by: 

AV(p) = ^TOoK 2 - UJ 2 2 ) P 2 Q{Z - Zq) 



(29) 



where Q(z) is the Heaviside function. The difference AV(p) is zero for the initial spherical 
configuration. 

Once AV(z) and AV(p) are diagonalized and added to the oscillator energy of the sphere 
+ ellipsoid system, which is 



E = huji[2n p + \m\ + v 1 + 1.5], 



(30) 



the level schemes of the two intersecting asymmetric oscillators with frequencies cu 1 and uj 2 
are obtained. 

C. Spin-orbit and orbit-orbit interactions 

The spin-orbit (1- s) and orbit-orbit (l 2 ) interaction terms generate the necessary single- 
particle level splitting to obtain the correct schemes of the individual fragments after sepa- 
ration. 
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The use of deformation dependent form of these operators has been introduced in [16| 
for the Nilsson model and in jTJJ] for the two center shell model; instead of the 1 operator 
one introduces: 



1 



VVxp 

m uj 2 



(31) 



where V is the asymmetric two-center oscillator potential. The usual expression for the two 
operators are: 



V(\-s) 



-2k%uj\ ■ s 



(32) 



Since one obtains the level schemes of two nuclei which lie in different mass regions, the 
strength parameters of the interactions k and fi will be different. The values we use for 
these parameters are: 



n n = 0.0588 n p = 0.0592 



(33) 



H n = 0.328 fi p = 0.335 
for actinide region, and for light nuclei region: 

n n = k p = 0.0601 
Mn = A* P = 0.448. 

Since the strength parameters are different for the asymmetric regions of the nuclear shape, 
they become ^-dependent operators as follows: 



(34) 



(35) 



K ■ hu(z) = Ki • TlUJi + (k 2 ' ^2 — K l ' hu>i)Q(z — Zq) 
KH ■ Huj(z) = Ki^i • TlUJi + (ft 2 /i2 ' ^2 — " hu>i)Q(z — Zq). 

To obtain a Hermitian operator for V(l ■ s) and V(\ 2 ) one has to use the anticommutator 



0: 



V(l ■ s) 
V(l 2 ) = 



«M-M^),(^F) : 



(36) 
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For the dependence of k and /x with respect to the elongation, we choose a linear dependence 
for the oscillator frequency along the z-axis: 

Ki = Kq+ — — (Kif - Kq) (37) 

Wif — UJq 

and the same law of variation for fi. Here i=l,2, k is the value for the initial nucleus and 
Kif for the final one. 

Finally one has to diagonalize the potential: 

AV(p, z) = AV(z) + AV(p) + Vi\ • s) + Vi\ 2 ) (38) 

together with the diagonal term of the two-center oscillator potential. The model provides 
the evolution from an initial level scheme toward two asymptotically independent single- 
particle schemes. With the introduction of large asymmetry between fragments, the shapes 
can simulate the existence of a soliton on the nuclear surface and assign it a microscopic 
structure. 



V. SHELL CORRECTIONS 

The level scheme of a soliton shape is used to obtain the shell corrections of the system. 
As the soliton is assimilated with an emerging fragment, it will provide the shell correction 
value of the independent nucleus of similar shape. Shell corrections are obtained by means 



of the Strutinsky procedure WR - One defines the shell correction energy as the difference 



between the total sum of the energy levels and a smoothed part of the spectrum: 

5E = J2 2E » - U. (39) 

V 

One calculates the smoothed part U with the help of a smoothed level density function g(e), 
which is obtained by averaging the real distribution g(e) over the whole energy spectrum: 

7^ I 7 / 
12 



where 7 = Y/huj and ((x) is the smoothing function. A common smoothing function is 
provided by 

1 m 

(( X )= e ~ x2 J2a2kH 2k (x), (41) 

where H n {x) are the Hermite polynomials. The coefficients a 2k are 

H 2k (0) 

a2k ~ 2*. (2*)r ( } 

The Fermi energy A of the smoothed level distribution is calculated as a solution of the 
particle number conservation: 



N p = 2 



[ X ~g(e)de. (43) 



Then, the total energy of the uniform level distribution U, reproducing the microscopic part 
which is not subjected to local fluctuations of the spectrum, is obtained as: 

r-A 



U = 2hw 



I 9(e)ede. (44) 
J — 00 



After performing the calculations, one obtains the following formula, which can be used 
directly 

5U = J2i^-erf(x Fv )} 

V 

H 7=r- • [2e u 22 a 2 kH 2 k-i + ia 2m H 2rn }}, (45) 



11 fc =i 

where erf(x) is the error function. Usually one chooses the upper order of the Hermite 
polynomials to be m=3. The variable xf v is given by 

x Fv = (46) 

7 

Shell corrections are calculated separately for protons and neutrons, and the results are 
added. 
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VI. RESULTS 



A first look at the energetic behavior of an emerging soliton is given in Fig. 2 in terms of 
the macroscopic energy surfaces. The LHS macroscopic energy surface corresponds to the 
formation of 40 Ca on the surface of 248 No, whereas the RHS represents the 16 0-like soliton 
on the surface of 224 Th. Variation along the elongation R corresponds to the increment in 
the soliton amplitude along the symmetry axis. A larger neck radius R% corresponds to a 
larger half-width L. The rear plane at R=0 is the spherical state of the system. Then the 
energy increases monotonously with a higher slope for small values of the neck radius. As 
i? 3 increases, the energy increase is smoothed by the necking. With the enhancement of the 
kinetic energy of the soliton the half-width becomes larger, except in the first stages of the 
process where the neck radius is very small. The ridge in energy has a maximum at the near 
touching spheres configuration for both reactions. The slope continues to increase for large 
R 3 , beyond the touching point value of the elongation R. 

The addition of shell corrections yields the total deformation energy shown in Fig. 3. As 
a first observation note the pronounced deformed ground state of 248 No as the first minimum 
in energy moves to R >0, and a much less but still deformed ground state for 224 Th. For 
both emerging solitons it is obvious that the energy path corresponds to large half-width 
values up to the top of the energy ridge; then they abruptly turn towards rupture point 
shapes (i?3=0). Hence, these potential energy surfaces suggest a three-dimensional curve 
as the path of minimum energy in the cluster emission. The potential barrier formed along 
the path of minimum energy values is obtained by minimization of the total energy in the 
multi-dimensional deformation space. 

The static paths, which a soliton with the internal structure of an emitted cluster has to 
follow, have been plotted on the contour maps of the energy surface in (R,R 3 ) coordinate 
space in Fig. 4. Again the LHS plot is the 40 Ca emission, and the RHS one corresponds to 
16 0. Apart from the first R values, where _R 3 is small, the solitons bypass the first energy 
peak by taking large neck radius values, i.e. large half-widths. As the energy increases on 
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the large R 3 side, the static path for both cases changes direction reaching the scission point 
where the clusters are emitted. 

The two barriers are plotted in Fig. 5 with a full line, together with the macroscopic 
energy (dotted line) and the shell corrections (dashed line) . One can see how the deformed 
ground state of 248 No is formed (LHS plot): due to shell corrections, the first minimum is 
at about 6.8 MeV of the total energy. This point becomes the ground state and the whole 
barrier in front of the emerging soliton is shifted with respect to this value. A two-humped 
barrier no higher than about 1.2 MeV blocks the 40 Ca emission. 

The situation is different for the emission of 16 from 224 Th. The ground state is only 
slightly deformed. A rather high one-hump barrier of about 11 MeV extends along the 
whole range of elongation R up to the scission point. Shell corrections decrease slightly the 
macroscopic energy values. The decrease is mainly due to the double-magic character of 
208 Pb which forms as the cluster emerges. 

One can state that 40 Ca-like solitons are energetically favored to form on the nuclear 
surface of a very heavy nucleus as 248 No. The formation of 16 0-like soliton on 224 Th is not 
energetically favored due to the high and large potential barrier it has to penetrate. 

The relative velocity distribution V of the two presumed solitons along the minimum 
energy path, together with the scaled values of the halfwidth L and the relative amplitude 
a = A/Ri, are plotted in Figs. 6. We investigated the evolution of these soliton parameters 
(as defined in section I), which are a function of the static energy evolution, parametrized 
by the distance between centers R. 

In the first stages, the tendency is that the amplitude and half-width increase with the 
elongation parameter, when the emitted cluster is emerging out from the parent nucleus 
(since their non-overlaping sector is increasing). During the formation of the cluster the 
half-width remains practically constant, since the surface energy controls this stage. When 
the two nuclei are well separated, the soliton envelope hardly fits the two spheres, and 
in this limit, the half- width approaches zero value. This gives the limiting configuration 
for this soliton model. These values of the half-width L (solid line) are compared with 
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those obtained analyticaly in || directly from the soliton amplitude, within the frame of 
the nonlinear liquid drop model (L-dashed line). We notice a good agreement for the half- 
widths within the range R ~ 4.5 — 14 Fm. The hydrodynamic soliton model is not valid 
anymore for separation parameter R smaller than 4-5 Fm, because of the dominating shell 
effects in this range. This can be noticed in a comparison between Fig. 2 and Fig. 3 for 
R < 4 — 5 Fm. For the first R values, the static paths follow the first energy peak, and 
jump from small toward large values for R%, Fig. 4, providing small half-widths (Figs. 6, 
L-solid line), while a pure hydrodynamic soliton would have larger half- widths for this range 
(L-dashed line). In the above range of validity of the soliton model, we calculate the relative 
velocity of the soliton (V^-dashed line), 0. The velocity is increasing with the amplitude of 
the soliton, hence with the elongation of the cluster-like emission shape. Fig. 6a displays 
the 40 Ca emission, and Fig. 6b represents the 16 emission. Soliton shapes at the begining 
and the end of the process are also shown. For lighter nuclei (like 16 0) the evolution of the 
parameters is smooth and monotonic. In the case of heavier nuclei ( 40 Ca) we obtained some 
oscillations in width and velocity, during the first half of the emission process, which can be 
related with the oscillations produced by the shell effects in the R3 parameter. 

This work was supported by the U.S. National Science Foundation through a regu- 
lar grant, No. 9970769, and a Cooperative Agreement, No. EPS-9720652, that includes 
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P. Draayer for a postdoctoral appointment in the Department of Physics and Astronomy, 
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Figure captions 

Fig.l Deformation space for two-center shell- model calculations. The neck radius R3 can 
vary from zero (intersecting spheres) to infinity (compact shapes). R (the distance between 
the two centers) and R 2 (radius of the emitted fragment) are also independent coordinates. 

Fig. 2 Macroscopic potential energy surfaces for 40 Ca emission from 248 No (left-hand-side 
plot) and for 16 emission from 224 Th (right-hand-side plot), as function of elongation R 
and neck radius R3. First maximum appears close to the touching point configuration in 
both cases (i? 3 =0). 

Fig. 3 Total potential energy surfaces (macroscopic plus shell corrections) for 40 Ca emis- 
sion from 248 No (left-hand-side plot) and 16 from 224 Th (right-hand-side plot). The de- 
formed ground state of 248 No is revealed as a minimum along R 3 axis at its origin. For 
both surfaces the closest energy maximum occurs at the tangent sphere configuration. The 
maximum energy value for larger i? 3 is not reached in the figure. 

Fig. 4 Contour plot of Fig. 3 with static path (dashed line) for 40 Ca emission (left- 
hand-side plot) and 16 emission (right-hand-side plot). As the elongation increases the 
amplitude of the soliton increases together with the half-width which is proportional to the 
neck radius R3. Once the touching point maximum is bypassed, both shapes decrease rather 
abruptly through necking towards scission. 

Fig. 5 The barriers along the static path for 40 Ca emission (left-hand-side plot) and 16 
(right- hand- side plot), together with macroscopic energy (dotted lines) and shell corrections 
(dashed lines). Shell corrections increase the total energy of the first energy minimum for 
248 No, thus the two- humped barrier is not higher than about 1.2 MeV. 16 emission from 
224 Th has a barrier of about 11 MeV. 

Figs. 6 The evolution of a = A/Ri, L (with shell corrections solid line, and without shell 
corrections dashed line) and V parameters in relative units, versus the elongation R in fm. 
The corresponding nuclear configurations are plotted for two situations: for the initial stage 
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when the emitted cluster is only slightly displaced off the common center, and the final stage 
when the two nuclei are almost separated. Fig. 6a displays the 40 Ca emission and Fig. 6b 
the 16 emission. In the later case, oscillations can be seen in the soliton parameters related 
to the shell corrections. 
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Fig. 2, Phys. Rev. C, Gherghescu et al. 




Fig. 3, Phys. Rev. C, Gherghescu et al. 
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